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Abstract 
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r*j ■ Two types of population models are well known - the continuous 

and the discrete types. The two have very different characteristics and 
methods of solutions and analysis. In this note, we point out that an 
iterative technique when applied to the continuous case mimics, sur- 
prisingly the discrete theory. The implication is that techniques and 
' conclusions of the latter theory can now be applied to the former case 

*0 ' (and vice versa). 
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. 1 Discrete Logistic Equation 

o 

Population growth in nature is seldom as smoothly continuous as a classical 
logistic curve suggests pQ. In a species with a short annual breeding season 
whose members live for several breeding seasons and die at any time of the 
year, a continuous record of population size would undoubtedly show sea- 
| sonal undulations. For many species in fact population growth is markedly 

discontinuous. These are, for example species whose members reproduce only 
once in their lifetime and die before their descendants' lives begin, that is 
each generation dies before the eggs are hatched, or the seeds are germinated, 
to form the successive generations. A continuous-time differential equation 
is then inappropriate as a representation of population growth. We have to 
consider instead difference equations. As this will be needed in the sequel, 
we touch upon the salient features. 
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1.1 First-Order Difference Equation 



One of the simplest systems an ecologist can study is a seasonal breeding 
population in which generations do not overlap. Many natural popula- 
tions,particularly among temperate zone insects (including the many eco- 
nomically important crop and orchard pests) are of this kind. In this situ- 
ation, the observed data will consist of information about the maximum or 
the average, or the total population in each generation [2J. The studies try to 
form a relation between the magnitude of the population in the generation n 
represented by x n and the magnitude of the population in generation n+1 
represented by x n+ i such a relation may be expressed in the general form 



The function F(x n ) will be what the biologist calls "density dependant", and 
a mathematician calls non-linear [3] . Equation (pQ) is thus a first order, non- 
linear difference equation. Equation (pQ) also describes many other examples 
in biology apart from population growth. 
For instance: 

1. In genetics, where the equation describes the change in gene frequency in 
time. 

2. In epidemiology where x n represents the fraction of population infected 
in time n. 

3. In economics where the relationships between commodity quantity and 
price are studied. 

4. In social sciences to study the propagation of rumors where x n could be 
the number of people to have heard the rumor after time t. 

In many of these contexts, and for biological populations in particular, there 
is a tendency for the variable x n to increase from one generation to the 
next when it is small, and for it to decrease when it is large. That is, the 
nonlinear function F(x n ) often has the following properties: F(0) = and 
F(x) increases monotonically as x increases through the range < x < 
A with F(x) attaining its maximum value at x = A, and F(x) decreases 
monotonically as x increases beyond x = A. A specific example is afforded 
by the equation 



x n+1 = F(x n ) 



(1) 




(2) 
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This is sometimes referred to as the logistic difference equation. Using 
the substitution x n = b/ax n Equation (j2J) can be written as 

F(x n ) = %n+l = dX n {\ Xn) (3) 

The behavior of the solutions of ([3]) is a function of the parameter 'a' Equation 
([3]) has meaningful solutions for 

< a < 4 

with x n measuring a non-negative quantity. 

Studies have shown [3] that the very simple nonlinear difference equation can 
possess an extraordinarily rich spectrum of dynamical behavior, from stable 
points, stable cycles to ultimately chaotic behavior. Thus the problem is far 
richer than the continuous case seen earlier. In this form equation ([3]) is a 
simple nonlinear difference equation. 



2 Examples of chaotic dynamical systems : 
The logistic map 

The discrete logistic map described by the single difference equation 

x n+1 = ax n (l - x n ) (4) 

as mentioned earlier, determines the future value of the variable x n +\ at 
time-step n+1 from the past value at time-step n. The time evolution of x n 
generated by this algebraic equation exhibits an extraordinary transforma- 
tion from order to chaos as the parameter a, which measures the strength of 
nonlinearity is increased [I]. 

Although nonlinear difference equations of this type have been studied ex- 
tensively as simple models for turbulence in fluids, they also arise naturally 
in the study of evolution of biological populations. 

For the purpose of illustration we consider the population of gypsy moths in 
the northern United States, which exhibits wild and unpredictable fluctua- 
tions from year to year. However we could equally well consider the evolution 
of economic prices determined by a nonlinear web model. 
Writing (J4]) in a slightly different form 

X n -\-l 0>X n &X, n 
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Figure 1 

Figure 1: Graphs of x n versus x n +i for different values of a 

we see that it is a simple quadratic equation, with the first term linear and 
the second term nonlinear. If the parameter a > 1, the population increases, 
if a < 1, the population decreases. If a > 1, the population will eventually 
grow to a large enough value for the nonlinear term — ax 2 nl to become 
important. Since this term is negative, it represents a nonlinear death rate 
which dominates when the population is too large. Biologically this nonlinear 
death rate could be due to the depletion of food supplies or the outbreak of 
diseases in an over-crowded environment. The dynamics of this map and the 
dependence on the parameter a which measures the rate of linear growth and 
the size of the nonlinear term, are best understood using graphical analysis. 
Consider the graphs of x n versus x n+ idisplayed in the Fig. ([I]) for four different 
values of a. Equation (@J defines an inverted parabola with intercepts at 

x n = and 1 

and a maximum value of 

Xn+i = a/4 
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at 

x n = 0.5 

Using these maps we can we can get a quantitative understanding of the 
dynamics of the logistic map in a quick way. 

Briefly, this graphical analysis tells us that if the normalized population starts 
out larger than 1, then it immediately goes negative , becoming extinct in 
one time-step. Moreover if a > 4, the peak of the parabola will exceed 1, 
which makes it possible for initial populations near 0.5 to become extinct in 
two time-steps. We will therefore restrict our analysis to values of a between 
and 4. 

For values of a < 1, the population always decreases to 0, as shown for 
a = 0.95 in Fig. ([I]) The intersection of the parabola with the 45° line at 
x n = represents a stable fixed point on the map. Because a is small a 
perturbation can be used to verify that almost all initial populations are 
attracted to this fixed point and become extinct. However for a > 1 this 
fixed point becomes unstable. Instead the parabola now intersects the 45° 
line at 

a — 1 

x = 

a 

which corresponds to a new fixed point. 

For values of a between 1 and 3 almost all initial populations evolve to 
this equilibrium population. Then, as a is increased between 3 and 4 , the 
dynamics change in remarkable ways. First the fixed point becomes unstable 
and the population evolves to a dynamic steady state in which it alternates 
between a large and small population. A time sequence converging to such 
a period-2 cycle is displayed in Fig. ([I]) for a= 3.2 : the population cycles be- 
tween two points on the parabola, x n ~ 0.5 and x n ~ 0.8, in alternate years. 
For somewhat larger values of a , this period-2 cycle becomes unstable and 
is replaced by a period-4 cycle in which the population alternates high-low, 
returning to its original value every four time-steps. As a is increased the 
longtime motion converges to period - 8,16,32,64, cycles, finally accumulat- 
ing to a cycle of infinite period a = a in f ~ 3.57 

Having observed a period doubling sequence in numerical experiments Feigen- 
baum was able to prove that the intervals over which a cycle is stable de- 
creases at a geometric rate of ~ 4.6692016. The tremendous significance 
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of this work is that this rate and other properties of the period-doubling 
bifurcation are universal in the sense that they appear in the dynamics of 
any system which can be approximately modelled by a nonlinear map with a 
quadratic extremum. Feigenbaum's theory has subsequently been confirmed 
by a wide variety of physical systems such as turbulent fluids, oscillating 
chemical reactions, nonlinear electrical circuits, and ring lasers. 
The investigation of period doubling in nonlinear dynamical systems pro- 
vides a superb example of the interplay between numerical experiments and 
analytical theory. However, this sequence of regular periodic orbits is only 
the precursor to chaos. Included below is a bifurcation diagram , showing 
the beginning of Chaos: 
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3 A Remarkable Mathematical Equivalence 



We would now like to deduce a mathematical equivalence between the con- 
tinuous and discrete cases and will comment on this. Let us consider the 
equation from the continuous growth 

dP 

— = f(P), setP^x (5) 

which gives 

x = [ f{x(t))dt (6) 



where the integral is over suitable limits. Here f generalizes the dependence 
of the right side of 

jv^Wj = r {M - P{t))P{t) r>0 (7) 

where P(t) is the population at given time t and M is the maximum sustain- 
able population [HIE] on the population P at that time. Let us solve equation 
(EJ) by the method of successive approximation ie., we try on the right side 
of ([6]) a tentative solution Xo(t). (EJ) can then be written in the form 

x x = F(x ) = J f(x )dt (8) 

where x± gives the next level of approximation. Before proceeding further, we 
remark that (jSJ) is in the form of the initial value problem,where the Lipschitz 
condition is necessary for the convergence of the iterative procedure [7]. In 
a similar manner we get from (jSJ) the more general equation 

x n+1 = F(x n ). (9) 

(Q can immediately identified with the discrete logistic equation for example 
(|3j). It must be stressed however that the discrete equation is based on a 
completely different foundation that is the subscript n in the discrete case 
represents the population at the n th generation, whereas in (Q x n represents 
the n th iteration or approximation of the population x = P of the continuous 
case. Nevertheless, this mathematical equivalence enables us to apply the 
conclusions of the discrete case including the domain of chaos. Thus chaotic 
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behaviour of x n of the discrete case would represent the lack of convergence 
of the iterates of the continuous case. In this specific example if 



F(x n ) = ax n (l - x n ), (10) 

then for a = 3.57 the above iterative procedure breaks down. In this case, 
using (flOj) . it follows from ([8]) that 

f(x)dt = ax(l — x). 

Finally we remark that for a more conventional approach to the above prob- 
lem reference can be made to Krempasky, [8]. 
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